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LAUNCHER  DYNAMICS  STUDY 


Part  M]^  PinsQVl 

(A  Two  Point  Mass  Moving  on  a  Flexible  Beam) 


ABSTRACT 

This  report  contains  an  approximate  analysis  of  the  motion  of  a  flexible 
launcher  rail  which  is  supported  by  a  flexible  understructure  and  Is  loaded 
by  a  missile  with  two  point  contact  which  moves  across  the  rail  under  the  action 
of  a  prescribed  thrust  force.  The  launcher  system  is  represented  by  a  uniform 
flexible  beam  which  is  supported  by  linear  springs  and  the  missile  is  constrained 
to  remain  in  contact  with  the  beam.  A  pair  of  coupled  integral  equations  are 
obtained  which  define  the  motion  of  the  missile  and  a  numerical  technique  of 
integration  is  developed  for  their  solution.  Finally  some  numerical  results  are 
presented  to  show  the  effects  of  some  of  the  parameters  on  the  motion  of  the  missile 


and  the  rail. 
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List  of  Symbols 


R  (t)  Thrust  force  on  the  missile 

W  Weight  of  the  missile 

m  Mass  of  the  missile 
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Other  symbols  are  defined  in  the  text  when  introduced. 


CHAPTER  1 


Development  of  Model 

A  typical  straight  rail  missile  launcher  mig^t  appear  as  sketched  in 
Pig.  1.  We  replace  this  with  the  model  shown  schematically  in  Fig.  2.  In 


F  igure  i 


the  undeformed  state  the  beam  is  straight,  uniform,  and  at  an  elevation  angle  6 
The  X  axis  is  taken  along  the  undeformed  axis  of  the  beam  and  the  y  axis  is 


'J 


Figure  2 


generally  downward.  The  beam  is  supported  at  two  points  by  linear  springs  and 
suitable  constraints  act  to  prevent  axial  motion  of  the  beam. 


2. 


We  assume  that  the  missile  is  a  rigid  mass  which  is  forced  to  move  along 
the  beam  under  the  action  of  the  prescribed  force  R  (t).  The  motion  of  the  mass 
along  the  beam  is  measured  by  the  coordinate  §  (t)  as  indicated  in  Fig.  2. 

During  the  motion,  two  points  of  the  missile  are  constrained  to  maintain  contact 
with  the  rail.  At  some  later  time  Tp  ,  with  §  (Tp  )  =  the  rear  shoe  also 
leaves  the  rail  and  the  missile  is  in  flight.  (It  is  possible  for  Tp  and  Tq  to  be 
equal).  We  refer  to  Tp  as  the  tipoff  time  and  assume  that  all  motion  occurs  in 
the  xy-plane  and  that  the  entire  system  is  initially  at  rest. 


CHAPTER  2 


Analytical  Develc^ment  of  the  General  Problem 
Consider  the  free  body  diagram  shown  in  Fig.  3. 


Figure  3 


Referring  to  Fig.  4,  let  9  be  the  rotation  of  the  missile  and  y  the 
displacement  of  its  center  of  mass.  The  equations  of  motion  of  the  mass  are 


4. 


If  we  let  Vp  and  be  the  transverse  displacements  of  the  points  of  contact 
then  for  small  displacements  and  rotations  of  the  missile 


L-L, 


y  -TT  yQ* 


and  therefore 


P=  m(g  h-JgVp  +  J^yQ) 

Q  -  m(gd  +  Jjyp-J2yQ) 

where 


L-L  j  L,<L-L,) 


(1) 

(2) 


J3  ' “Vz  ^ r=  geos  eg 


In  order  to  obuin  a  good  approximation  of  (  (t)  we  integrate  the  equation 


5. 


mg  =  R  :  g(0)  --0.  g(0)=  g^^ 

where  R  is  assumed  to  be  a  known  function. 

For  the  beam  of  Fig.  3,  the  equation  of  motion  is 


gj  +  pA  =  P6(g,x)  + Q6(g  +  L,x)  (3) 

8  X  8  t 

where  6  is  the  Dirac  Delta  or  Unit  Impulse  function  defined  in  the  usual 
manner,  i.  e. 


X 


1 


/ 


6  (g,x)  dx  = 


0,  x^  <X2<S 

i  1,  Xj  <  g  -c  x^ 
0,  g  <:  Xj^  ^  x^ 


Since  the  system  is  initially  at  rest,  the  rail's  Initial  velocity  y^  (x,  0)  is  zero 
and  its  initial  displacement  y  (x,  0)  is  due  to  the  static  loads  produced  by  the 
mass  m  . 

Assume  that  the  solution  of  Eq.  (3)  may  be  written  in  the  form 


y(x,t) 


(«) 


n=J 


-th 


where  (x)  is  the  n  normal  mode  for  the  spring  supported  beam  of  Fig.  3, 
and  each  q^  (t)  may  be  regarded  as  a  generalized  coordinate.  The  set  of  functions 


are  the  eigenfunctions  of  the  equation 


d*X^(x) 


-X„X„(x)  =  0. 


n  E  I  n 


n  =  1,  2,  3, . . . 


(5) 


6. 


subject  to  the  appropriate  bounda:  y  conditions  for  the  beam.  (The  detailed 
solution  of  this  equation  is  given  in  Aj^ndix  A),  is  the  n^^  natural  frequency 
for  the  beam.  The  set  of  functions  (x^  also  have  the  property 


J 


(X)  X„  (X)  dx  = 
m  n '  ' 


{0,  m  /  n 
1,  m  =  n 


(6) 


Having  determined  the  functions  Xj^  (x)  we  can  e}q)and  the  Dirac 
functions  in  the  following  series 


eo 

M§.x)  =  ^X^(g)Xj^(x) 
oo 

6(g+  L.x)»  2Xj^(e+  L)Xj^(x) 

n?l 

Using  Eqs.  (4^,  (7),  and  (8),  Eq.  can  be  written 

CD 


(7) 

(8) 


nn  sn 


and  with  the  use  of  Eq.  (5)  this  e3q>ression  becomes 


(X) 


Upon  comparing  coefficients  of  the  functions  ^X^  (}o|  ,  we  obtain 

■  pT  ■  n- 1,2.3,... 


0<t<Tf. 


(9) 


Since 


y(x.0)=  2, 

y.  (x.  0)  =  1  XU)L(0)-0 


7. 


it  follows  that  the  solutions  of  Eq.  (9)  are  subject  to  the  conditions 

q„  (0)  =  f  X  (x)  y  (X,  0)  dx  (10) 


\(0)  =0 


(11) 


(An  e:q>ression  for  q^  (0)  is  derived  in  Appendix  A).  Thus  the  solution  of 
Eq.  (9)  becomes 

(t)  =  q^  (0)  cos  +  [g  (t)|  +  Q  (T)  [g  (t)  +  ^jsin  Ujj 


+  Q(t)X, 

n  ~  1|2|3. 


(t-T)  dr 
(12) 


Upon  setting  x  =  g  (t)  in  Eq.  (4),  we  obtain  the  deflection  of  the  point 
of  contact  of  the  rear  shoe: 


n=l 

Similarly  for  x  =  g  (t)  +  L,  we  find 

Yq  (t)  =  y(g  +  L,t)  =  5  Xj^  [g  (t)  +  l]  q^(t) 


(13) 


(14) 


n=] 


Then  combining  Eqs.  (1),  (2),  and  (14)  and  substituting  this  result  into  Eqs.  (1^ 
and  (1^  we  obtain  the  pair  of  coupled  integro-differential  equations  in  the 

shoe  deflections  and  their  accelerations,  i.  e. 

.t 


yp  (t)  =  F  ( t )  +  J* jkj  (t, 7)  +  K2  (t, t)  y'p  (t)  +  Kj  (t, T)  jTq  (T)j 


d  T 


0  <  t  <  X 


(15) 


(t)  =  C{t)+J jic^  (t,  t)  +  Kj  (t,  t)  Vp  (t)  +  (t,  t)  Vq  (T)j 


Q 

d  T 


0<  t<T,r 


(16) 


8. 


where 

F  (t)  =  [e  (t^  (Q)  cos  t 

G  (t)  =  ^  Xjj  [g  (t)  +  l]  (0)  cos  t 


K2(t,T)=  p 


K, 


m 


|hX„[UT)]tdX„[S(T)  +  L]j 


Kg  (t,  T)  = 


{jl  *n  [«")  +  ■■]  -  J3  Jtn  N} 

Note  that  the  (t,  r)  are  non- symmetric  kernels  with  the  property  (t,t )  =  0  . 

In  order  to  put  Eqs.  (15)  and  (16)  in  a  more  tractable  form,  let 


Up  (t)-  Vp  (t) 

UQ(t)*'yQ(t) 


Then  since  the  missile  is  initially  at  rest, 
Yp  (t)  =  ^ [t  -  7)  Up  (t)  d  T  +  yp  (0) 
Yq  («)  =  At  -  T)  (T)  d  T  +  Yq  (0) 

oJ 


(21) 

(22) 


9. 


in  which  Vp  (0)  and  (0)  are  the  initial  displacements  of  the  shoes  and  are 
con4>uted  from 


(0)=  |^q„(0)X„(f„) 


yq  <0)  = 


Therefore  Eqs.  (15)  and  (16)  can  be  written  in  the  form 

Yp  (0)  =  F  (t)^+ 1* (t,  t)  +  (k2  (t,  T)  -  (t-T)]  Up  (T)  +  Kg  (t.  T)  Uq  (T)j  d  T  (23) 
Yq  (0)  =  G  (t)  t)  +  Kg  (t,  t)  Up  (t)  +  [Kg  (t,  t)  -  (t-i)J  Uq  (t)^  d  t  (24) 

These  equations  are  to  be  solved  for  Up  (t)  and  (t).  (In  the  next  section  we 
will  present  a  numerical  technique  for  obtaining  appropriate  solutions).  The 
motion  of  the  missile  is  then  computed  from  E  s.  (21)  and  (22)  and  the  motion 
of  the  beam  can  be  obtained  from  Eqs.  (4)  after  P  and  Q  are  determined  from 
Eqs.  (1)  and  (2)  and  qQ(f)  Is  found  from  Eq.  (1^. 

Eqs.  (23)  and  (2^  are  applicable  when  both  shoes  of  the  missile  are  in 
contact  with  the  rail.  When  the  front  shoe  loses  contact  with  the  rail,  Q  becomes 
and  remains  zero.  For  this  case  we  set 


y  (x.  t  )  = 
obtaining 


tin 


(t) 


Tq<  t<Tp 


(25) 


JT  Xn«) 


n  =  1,2,  3, ... 


(26) 


10. 


subject  to  the  continuity  conditions  (Tq  )  =  <1^  (Tq  )  and  (Tq)  =  q^^  (Tq  ), 
1.  e.  the  displacements  and  velocities  of  the  beam  are  continuous.  Therefore, 

q*n  (*)  =  <'^Q  >  “n  ^*n  <  V  ‘"n  > 


■K/ 


sin  (t  -  t)  d  T 


(27) 


n  *  1,  2,  3, .  •  a 

TQ«it<Tp 


where 


P*  =  m 


Jl  +  hd 


[i-?J 


(28) 


is  found  by  setting  Q  =  0  in  Eq.  (2)  ,  solving  for  y^  ,  and  then  substituting 
this  result  into  Eq.  (1). 

Differentiating  Eq.  (1^  twice  we  obtain 

Vp  <«  '  Z[*n  <*>  '>"<■>*  ^*0  «>  %  <•»] 


n=l 


0  <  t  <  Tp 


and  similarly,  from  Eq.  (25)  after  setting  x  =  $  (t) 


^’“P  <■>  '  £Pn  «>  ■>'»  +  2  X„  (5)  (0  .  X„  (t)  (t)] 


Tq<«<Tp 


If  we  set  t  =  ’^Q  these  equations  and  subtract,  we 


find 


where  use  has  been  made  of  Eqs.  (1),  (2),  (9),  (26),  and  (28)  to  simplify  the 
result. 

In  Eq.  (25)  we  now  substitute  g  (t)  for  x  and  proceed  as  before  to  obtain 
the  displacement  of  the  rear  shoe, 


y*p  (t)  =  F* 


(t)  +  J[K*i(t,Tj][g-'y*p(l)J  dT 


(30) 


Tq<  t<Tp 


where 


F*(0  = 


sin^„ 

I  I  n 


(t,T)  =  m 


CO 

J^  +  hd  ^ 


V''‘*n<V  cosu.^(t-TQ) 
’‘n  [«  (^>] 


^  X 


n  [« 


<0 


sinw^(t  -  t) 


(K*^  (t,  t)  is  also  a  non-symmetiic  kernel  with  the  property  K*^  (t,t)  =  0  ). 
If  we  now  set 


u*p(t)  .  y*p(t) 


then  upon  integrating. 


y*P  ^  pP  'Z  (*  ■  Tq)  ^  yp  (Tq)  -  J  (t*T)  [g*u*p  (t^  dr 


Therefore  Eq.  (30)  can  be  written  in  the  form 


('  -  V  [yp(TQ)  t  ‘  (t  -  Tq)  *]  t  yp  (Tq)  =  F*  (t)^^ [k"j  (t-,^[g-u*p(T)J 


dr 


Tq<  t  <Tp 


(31) 


12. 


For  the  condition  of  only  the  rear  shoe  in  contact  with  the  rail  we  obtain 
a  single  integral  equation  which  we  must  solve  for  g  -  u*p  (t).  (In  the  next 
section  we  will  present  a  numberical  technique  for  obtaining  appropriate  solutions). 
The  motion  of  the  system  is  now  computed  as  before  except  that  the  determination 
of  is  replaced  by  the  condition  Q  =  0. 


CHAPTER  3 


13. 


Numerical  Solution  of  Equations 

An  approximate  solution  of  the  equations  of  motion  can  be  obtained  in 
the  following  manner.  For  definiteness,  we  suppose  the  thrust  curve  is  as  shown 
in  Fig.  5. 

R.(t) 

\4AX 


Figure  5 

If  t^  is  the  thrust  buildup  time  and  fhe  value  of  R  (t)  at  t^  ,  then 


Now  we  divide  the  interval  0  <  t^  into  M  equal  intervals  of  time 
T 

At  :  in  length.  Observing  that  Up(0)  =  (0)  =  0,  (t^t )  -  0,  and  applying 


14. 


the  trapezoidal  rule  of  Integration,  we  obtain  the  following  approximations  for 
Eqs.  (15)  and  (If). 


y  (0)  -  F  (N  At  )  ^ 

^  i  (N  At  .  0  )  2  <  N  <  M 


At 
N-2 


+  y  Gc,(N  At.  j  At)+  |kj(N  At.  j  At)  -(N-J)  AtJ  Up(i  At) 

+  K3(N  At,  j  At)  (j  At)V 

[n  At,  (N-1)  A^  +  .^2  ^*1  Up^N-1)  At] 

+  Kg  [n  At,  (N-1)  At]  Uq  [(N-1)  At] 


(32) 


yQ(0)-G(NAt)  ^ 


At 

N-2. 


I  K^(NAt,  0)  2  <N 


.M 


+  ^^4  j  At,  j  At  )  Up  (j  At) 

+  |kg(NAt,  j  At)  -(N  -  j)  At]  UqO  At)| 
+  [n  At,  (N-1)  At]  +  Kg  [n  At,  (N-1)  At]  Up  [(N-1)  At] 

+ 1*^6  ^*1 


(33) 


We  begin  the  numerical  process  with  N  =  2  and  then  use  successive  substitutions 
to  carry  out  the  numerical  integration  to  N  =  M.  The  advantage  of  using  this 
form  for  the  equations  is  that  they  are  recursive  in  nature  and  therefore  do  not 
involve  the  Inversion  of  a  matrix  for  their  solution.  The  bulk  of  the  computations 
will  be  in  computing  the  kernels  K^  (t,  r)  at  the  various  sub- intervals  of  time. 


15. 


The  sums  on  j  create  no  large  amount  of  computation  because  the  kernels  can 
be  rewritten  so  that  the  corresponding  values  at  J  +  1  are  obtained  from  those 
at  j  by  adding  a  term. 

Because  (t,  t)  =  0,  we  use  Eqs.  (32)  and (33)  to  compute  Up(t)  and 
lyj  (t)  only  up  to  Tq  -  At.  In  order  to  compute  the  values  at  Tq  a  second  degree 
polynomial  is  fitted  to  the  preceeding  three  values  resulting  in  the  formulae 


Up(Tq)  'Up  [(M-3) 

Uq  (Tq)  =  Uq  [(M-3) 


^t]  +  3  |up  [(M-1)  At]  -  Up  [(M-2)  At]| 

A  ]  +  3|uq  [(M-1)  At]  -  Uq  [(M-2)  a  t]J 


The  main  factor  influencing  the  choice  of  the  number  of  subdivisions 

To 

M  =  is  the  factor  sin  (t  -  t)  =  sin  (N-j)  At  in  each  kernel.  In 

numerical  computations  we  have  to  truncate  the  series  expressions  at  some 
value  n  =  p.  Then  the  factor  sinUp  (t  -  t)  could  range  over  many  complete 
cycles  in  the  interval  0  <  t  <  Tq  .  Wp  (n-j)  A  t  changes  by  the  amount  Wp  At 
with  each  increase  in  j  and  prudence  demands  that,  at  the  very  least,  Wp  At  < 
(In  the  computations  made,  the  worst  case  had  five  points  per  cycle  in  the  third 
mode  ,  i.  e.  p  =  3  was  usei). 

For  Eq.  (31)  we  divide  the  interval  T^..^  t  ^  Tp  into  M*  equal  intervals 

'P  -  'P  ^  * 

of  time  At*  =  ^^4,^  in  length.  Observing  that  K*|^  (t,  t)  =  0  ,  the 

trapezoidal  rule  of  integration  applied  to  Eq.  (31)  yields 


lOl  =1 


16. 


N*  At*  [yp(TQ)  +  ^N*  At*  g]  +  yp(TQ)  *  F*  (Tq  +  N*  A  t*  ) 


2<N*^  M* 


=  ^  [k*i  (Tq  +  N*  At*.  Tq)  +  N*  [g  -  u*p  (Tq  )]  (34) 

N*-2 

+  2  [k*^  (Tq  +  N*  a  t*.  Tq  +  j  At*)  +  (N*-J)  At*]  [g-u*p  (Tq  +  J  At*)] 

j=l 

+  |k*^  [Tq  +  N*  At,  Tq  +  (N*-1)  a  t]  +  At*||g  -  u*p  [Tq  +  (N*-l)  At]j 


where  we  have  set  t  =  Tq  +  N*  A  t*.  We  begin  the  numerical  process  with 
N*  =  2  and  carry  out  the  same  resubstitution  process  that  was  used  in  Eqs.  (32) 
and  (3^ .  Since  Eq.  (34)  can  be  used  only  up  to  Tp-  At*,  we  again  compute 
the  value  at  Tp  from  a  second  degree  polynomial  approximation.  The  time 
interval  At*  is  chosen  comparable  to  At  by  taking  M*  as  the  next  integer  larger 
than(Tp-TQ)/MTQ  . 

Having  determined  the  values  of  Up,  Uq  ,  and  g  -  u*p  at  each  of  the 
sub-interval  points,  we  are  in  a  position  to  determine  values  for  the  following 
parameters.  Each  of  the  indicated  integrations  may  be  carried  out  using 
Sin^son's  rule  for  example,  and  only  p  normal  modes  are  considered  in  the 
series  e:q)ansi(m. 

For  ^  ,  the  pitch  angle  of  the  missile 


^(t)  = 


ft  (0)  +  J  (t-T)  (t)  -  Up  (T)]  d  T  ,  0  ^  t  <  Tq 

0*'  M 

[(t  -  Tq)  S'  (Tq)  +  ?  (Tq)  +  ,-X  /  <‘'^)  [«  -  “’p  <^)]  ^ 

2  m  W 


t<Tp 


For  ^  ,  the  pitch  velocity  of  the  missile 


Vi  !• 
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•  [i J[“q ■  "p «]  " ’• 

0  (t)=  <  (r  ^ 

■  /  [i  -  u*p  (t)] 


d  T 


,  0^  t<T^ 


t<T„ 


For  y  ,  the  center  of  gravity  displacement  of  the  missile 


dT,  0<t<T^ 


y(t) 


p(0)  +  [h  Up  (t)  +  d  Uq  (t)] 

j(t-TQ)f (TQ)  +  y(TQ)t^(t-TQ)^i-  r (t-T)  [g-u*p(T)] 


Tq«  t<Tp 

For  y  ,  the  center  of  gravity  transverse  velocity  of  the  missile 


(t) 


I  [h  Up  (t)  +  d  u  (t)1  d  t  ,  0  ^  t 

f  Tr-  1 

V  .  ^  V  -  1  J  L«  -  “*p  " 


|y(TQ)  +  (t-TQ)g-— ^ 


Q 

dT,  Tq<  t<Tp 


For  yjjp  ,  the  tip  displacement  of  the  rail 

P 


(c)  (t) 

|jx„(c)q*„(0 

For  y,pjp  ,  the  tip  transverse  velocity  of  the  rail 


,  0  <  t  < 


''tip  ' 


Tq<  t<  Tp 


PripW'  < 


I; 


(t) 


O^t^Tr 


Tq^t^Tp 


For  y' jjp  >  the  tip  angular  displacement  of  the  rail 


d  T  , 


.  0<t<TQ 
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X'„(c)qn(t) 


y  TIP  ~  ' 


n=l 


X;(c)q„(.) 


Tq<  t<  Tp 


^  TIP  ’  ®“gular  velocity  of  the  beam 


y’-np  <‘)  = 


0  <  t  <  T. 


t<T, 


CHAPTER  4 
Numerical  Results 
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Table  1  gives  the  nominal  values  of  the  parameters  specifying  a 
particular  system. 

Table  1 


Parameter 

Description 

Value 

Rear  spring  support 

CO 

^2 

Front  spring  support 

9,  230  Ib/in 

a 

X  coordinate  of 

46  in 

b 

X  coordinate  of  k2 

160  in 

c 

Total  rail  length 

328  in 

El 

Flexural  rigidity  of  rail 

54  X  10^  Ib-in^ 

pA 

Mass  per  unit  length  of  rail 

0.018325  Ib-secVin-in 

®E 

Angle  of  elevation  of  rail 

0.  8  rad 

L 

Length  between  shoes  of  missile 

144  in 

Length  from  rear  shoe  to  center  of 

gravity  of  the  missile 

100  in 

W 

weight  of  missile 

4150  lb 

J 

Prlnc4>al  transverse  moment  of 

inertia  of  missile 

41,  200  lb*sec^in^/in 

^o 

X  coordinate  of  rear  shoe  at  t  »  0 

4  in 

*^MAX 

Maximum  thrust  force  on  missile 

124,  686.  6  lb 

‘‘l 

Length  of  guidance  for  two  shoe  contact  180  in 

Total  length  of  guidance 

180  in 

*0 

Thrust  build-up  time 

0. 100  sec 
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The  following  variations  of  the  nominal  value  s  were  also  considered 

Table  2 


System 
Soft  beam 
Stiff  beam 
Soft  spring 
Stiff  spring 
Short  overhang 

Short  tipoff 
Long  tipoff 
Lig^t  beam 
Heavy  beam 


Parameters  Changed 

®  ■  ■  -5  )nom 
=  I  =  -5  nom 

^2  ""  2  ^^2  ^  nom 

^2  ""  ^^2  ^  nom 

b  =  (c  ) 

'  '  nom 

^2  ”  ^  c  ^  nom  ^^2^  nom 
‘*2  =  '-0'<‘'l>nom 

pA4<PA)„om 

"*=  2(pA)„„„ 


For  the  nominal  system  results  were  computed  using  one,  two  and 
three  modes.  The  results  for  and  ?  are  shown  in  Figs  6  and  7.  For  the 
parameters  involved,  two  modes  apparently  furnish  values  as  accurate  as  three. 

Figs  8  through  15  present  graphically  the  results  obtained  for  the  systems 
given  in  Tables  1  and  2.  Three  modes  were  used  in  all  of  the  computations. 

The  effect  of  long  and  short  tipoff  are  given  in  Table  3. 
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Table  3 


Nominal  Short  Tipoff  (Tp  = 

222.7) 

Long  Tipoff  (Tp 

e 

3.300 

10. 320 

« 

24.960 

199.  210 

y 

0.061 

0. 142 

• 

y 

0. 944 

1.742 

y-np 

0. 074 

.057 

• 

y-rip 

0.342 

-0.513 

y'Tip 

3.  340 

2.670 

y'xip 

4.  690 

2.  770 

The  natural  frequencies  of  the  unloaded  beams  in  the  above  systems  are 

also  of  interest. 

These  are  presented  for  the  first  three  modes  in  Table  4. 

Table  4 

System 

Natural  Frequency,  cycles/sec 

fl 

^2 

^3 

Nominal 

4.6 

52.6 

156.3 

Soft  beam 

4.6 

37.7 

110.5 

Stiff  beam 

4.7 

73.8 

221.0 

Soft  spring 

3.3 

52.2 

156.3 

Stiff  spring 

13.0 

58.9 

156.  3 

Short  overhang 

4.7 

52.1 

156.3 

Short  tipoff 

4.6 

52.6 

156.3 

Long  t4>off 

4.6 

52.6 

156.3 

Li^t  beam 

6.6 

74.3 

221.0 

Heavy  beam 

3.3 

37.2 

110.5 

283. 1) 
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t,  Milliseconds 


Mils 


Figure  8 
^  vs.  t 


4.  57 


t,  Milliseconds 


Mils/Second 
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t.  Milliseconds 


Feet 
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,  Feev  Second 


t,  Millifeconds 


,  Mils/Secxjnd 
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t,  Milliseconds 


32. 


APPENDIX  A 
Solution  for  Eq. 

Fig.  A.  1  shows  the  beam  of  Fig.  3  in  an  unloaded  state. 


y  F  igure  A.  1 

th 

The  n  normal  mode  (x)  for  this  spring  supported  beam  are  the 


eigenfunctions  of  the  equation 


whose  general  solution  is 


(A.  1) 


X„  (jO  =  A'  sin  X  +  B'  cos  x  +  C’  sinh  x  +  D’  cosh  \  x  (A.  2) 
n  n  n  n  n  n  n  n 

Since  both  ends  of  the  beam  are  free,  we  must  have 


X"n  (0)  -  X'"^  (0)  =  X^ •’  (c)  =  (c)  =  0  (A.  3) 

At  the  supports,  i.  e.  at  x  =  a  and  x  =  b,  X^  (x),  X'^  (x)  ,  and  X"^  (x) 
must  be  continuous  and 
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E  I  X'”jj  (a‘ )  -  E  I  X"’jj  (a^  )  ’  *^1  ° 

E  I  X"'^  (b‘  )  -El  Xjj  "  (b"'" )  -  k2  \  ° 

In  view  of  Eqs.  (A.  4)  and  (A.  5)  we  make  the  following  definitions 
0<x<a:  X^  (x)  =  (^0  =  sin  (a-x)  +  cos  (a-x)  +  sinh  (a-x) 

+  Dj^  cosh  Xjj  (a-x)  (A.  6) 

a<x<b:  Xjj  (x)  =  (x)  =  E^^  sin  Xj^  x  +  cos  X^^  x  +  sinh  X^^  x 

+  H  cosh  X  X  (A.  7) 

n  n 

b<x<c:  Xj^  (x)  =  Tij^  (:0  =  I^  sin  X^^  (x-b)  +  cos  X^^  (x-b)  +  sinh  X^^  (x-b) 

+  L  cosh  X„  (x-b)  (A.  8) 

n  n 

Satisfying  the  conditions  on  X^^  (30  at  the  supports  and  results  in 

A„  =  -E„  cos  X„  a  +  sin  X„  a  +  o  A  u  (a) 
n  n  n  n  n  n  n 

B  =  E  sin  X  a  +  F  cos  X„  a 

n  n  n  n  n 

C  =  -G  cosh  X  a  -  H  sinh  X  a  -  a  A„  n  (a) 
n  n  n  n  n  n  n 

D„  *  G„  sinh  X„  a  +  H  cosh  X„  a 
n  n  n  n  n 

In  =E^cosX^b-F^8inX^b-^A^M„(b)  (A.9) 

J  =  E  sin  X  b  +  F  cos  X  b 

■’n  n  n  n  n 

K  =  G  cosh  X  b  +  H  sinh  X„  b  -  /3  A„  n  (b) 

11  n  n  n  n  n  n 

L  =  G  sinh  X  b  -i-  H  cosh  X  b 
n  n  n  n  n 
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where 

k  k 

1  *1  .  .  1  *2 

2  "El  P  ■  I  'FT 

Using  Eqs.  (A.  9)  we  now  have 

0<x<a:  (x)  =  <|>^(x)  =  ji^(x)  +  a  ^  H^(a)  |8ln  {a-x)  -  slnh  X^^  (a-j^  .  (A.  10) 

a^x^b:  X„  (x)  =  u,  (x)  =  E„  sin  X„  x  +  F„  cos  X^  x  +  G  slnh  X  x  +  H  cosh  X„  x  (A.  11) 
^  ^  n  nnnnn  nn  n  ' 

b^x^c:  (x)  =  ii^(x)  =  n^(x)  -  /3  [sin  X^^  (x-b)  +  SlnhXj^  (x-b)]  (A.  12) 

From  the  end  conditions  of  Eq.  (A.  3)  we  then  have 


E  a  A  sin  X  a  (sin  X  a  +  slnh  X_  a  ) 
n  n  n  n  n  ' 

+  F  ll  +  a  A  cos  X  a  (sin  X  a  +  slnh  X„  a  )1 
nL  n  n  n  nu 

+  G.  aA  slnhX  a  (sin  X  a  +  slnh  X  a  ) 
n  n  n  '  n  n 

-t-  H  fa  A  cosh  X^  a  (sin  X  a  +  sinh  X  a  )  -  1 1  =  0 

nLn  n  n  n'j 


E„  Fa  A„  sin  X„  a  (cos  X_  a  +  cosh  X„  a  )  -  1 1 
nLn  n'  n  n  J 

+  F„  «  A„  cos  X„  a  (cos  X„  a  +  cosh  X„  a  ) 

Q  n  n  n  n 

+  Fl  +  a  A„  slnh  \  a  (cos  X„  a  +  cosh  X  a  )1 
nL  n  n  n  n'j 

+  a  Ajj  cosh  Xj^  a  (cos  X^  a  +  cosh  Xjj  a  )  =  0 
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sin  c  +  /3  Ajj  sin  b  jsin  X^^  (c-b)  -  sinh  X^^  (c-b)|| 

+  Fj^  ^-cos  c  +  /3  Ajj  cos  X^^  b  jjsin  X^^  (c-b)  -  sinh  X^^  (c-bj|| 

-K;  ^inh  X  c  +  p  A  sinh  X  b  fsin  X  (c-b)  -  sinh  X  (c-b)]/ 
n\,n  n  nLn  n.U 

+  ^cosh  c  +  /3  Ajj  cosh  X^^  b  ^sin  X^^  (c-b)  -  sinh  X^^  (c-b)j| 

®n  ^  ^n  ^  ^  ^  ^n  \  ^  If°®  ^n 

+  F  <sin  X  c  +  p  A  cos  X„  b  Icos  X„  (c-b)  -  cosh  X  (c-b)!^ 

n\^  n  n  nL  n  -JJ 

Gn  ^cosh  \^c+  p  A^  sinh  X^^  b  jjcos  X^^  (c-b)  -  cosh  X^^  (c-b)j| 

H  <sinh  X  c  +  /3  A  cosh  X  b  fcos  X  (c-b)  -  cosh  X  (c-™/ 

n\  n  n  nLn''  n  jJ 


=  0 


=  0 


These  last  four  equations  are  a  homogeneous  set  and  for  non- zero  values  of 

E  ,  F  ,  G  ,  and  H  the  determinant  of  the  coefficients  must  be  zero.  The 
n  n  n  n 

roots  Xq  of  this  determinant  are  the  eigenvalues  for  the  problem  . 

We  solve  the  homogeneous  equations  in  the  form 


E„  =  e_  z^ 
n  n  n 


F  =  f  2 
n  n  n 


^n  '  ^  *n 


Hn  =  K  K 

n  n  n 


where  at  least  one  of  the  e  •  f  .  g_  ,  and  h  is  one. 

n  n  ”p 


of  z  such  that 
n 


Then  we  dioose  the  value 
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/x„wx„(x)dx=  ‘  [xx„^  ’  x„x;-  ^ 

-  X'  X"’T  =  1 

For  this  problem  this  means  that  is  chosen  such  that 

\  ^n  -4aari'^(a)Ti^(a) 

+  6  /3  [%  (b)]  -  4  ^  b  ^’jj(b)  n^(b)  +  c  (c^  ^ 


=  1 


To  prove  that 

'C 


J 


X  (x)  X  (x)  dx  =  0  for  m  /  n 
m  n 


we  write 


X  x^v 

-  X  ^X  X 

=  0 

m  n 

n 

m  n 

X 

X 

X 

1 

=  0 

n  m 

m 

m  n 

Subtracting  and  integrating  over  the  length  of  the  beam  results  in 


/*■ 


X^JJ.  )  d  x 


Integrating  the  right  hand  side  by  parts  yields 
>c 


>/ *m  «*„««*■  [♦„  *"X  *  t"”- ""  J. 

^  c  a 

^  ”  m  ”■•  Jb  -  [♦.♦m Jo 
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And  upon  applying  the  conditions  on  (x),  (x)  ,  and  we  have 

'C 


(X  jj  -  X.  m 


J 


4  )  I  (X)  X„  (X)  d  X  =  0 


m  '  n 


4  4/ 

Since  in  general  ^  ~  ^  xn  '  ^ 


/-■ 


(3^  Xjjj  (x)  d  X  =  0  for  m  /  n 


Thus  proving  the  orthogonality  of  the  beam  functions 
For  Eq.  (10)  we  write 


{x„(xj}  . 


0 


(x)  y  (x,  0)  d  X 
Integrating  by  parts  results  in 


JXj^(x)  y*^(x.( 


+ 

0 


,0)  dx 


IV 


For  a  beam  statically  loaded  with  concentrated  loads  y  (x,  0)  is  zero  and 
for  the  beam  of  Fig.  A.  1  y  (x,  0),  y'  (x,  0)  ,  and  y"  (x,  0)  must  be  continuous 
and  y'"  (x,  0)  satisfies  the  same  discontinuity  conditions  which  X^(30  satisfies. 
Therefore  the  final  results  for  (0)  can  be  shown  to  be 


,^(0).  2!^  ).dX,<|  .L)] 

P  A 


(A.  13) 
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